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The gap in semiconductor nanocrystaUites has been extensively studied both theoretically and 
experimentally over the last two decades. We have compared a recent "state-of-the-art" theoretical 
calculation with a recent "state-of-the-art" experimental observation of the gap in Si nanocrystallite. 
We find that the two are in substantial disagreement, with the disagreement being more pronounced 
at smaller sizes. Theoretical calculations appear to over-estimate the gap. Recognizing that the 
experimental observations are for a distribution of crystallite sizes, we proffer a phenomenological 
model to reconcile the theory with the experiment. We suggest that similar considerations must 
dictate comparisons between the theory and experiment vis-a-vis other properties such as radiative 
p^ . rate, decay constant, absorption coefficient, etc. 
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I. INTRODUCTION 



Semiconductor nanocrystallites, more popularly known as quantum dots (QDs), have been extensively studied over 
the past decade and a half. The system is interesting from the point of view of basic physics, with the carriers 
being confined to an essentially "zero dimensional" structure. The efficient luminescence observed in some of these 
crystallites makes them promising candidates for opto-electronic devices. Further, the inexorable drive towards device 
r»^ ,: miniaturization makes them technologically significant. i-.|-. 

f^ ■ The earliest theoretical works in this field were reported in early 1980'sEl'H. For over a decade after, most of the 

^-H [ theoretical works reported were based on effective mass theory (EMT) and tight-binding semi-empirical approaches. 
^^ I These works predicted the experimentally observed trends for quantum confinement, i.e. the inverse dependence of 
enhanced band gap on the nanocrystallite size. These calculations were, however, performed for a single nanocrys- 
tallite, whereas experiments are performed on an ensemble of QDs of varying sizes. We have earlier pointed out 
Ch ' the existence of an ensemble of QDs of varying sizes must be taken into account in any theoretical formulatio: 
rj—t , Though improved theoretical (calculations were pursued later, this aspect has largely been ignored. 
^ • In a recent work Ogut et alB reported a "state-of-the-art" theoretical calculations based on first principles. This 

Q ' work is claimed to be in excellent agreement with an early experimental workQ and is claimed to be superior to semi- 
O • empirical calculations. We compared their calculations to a later "state-of-the-art" experimental results of Buuren 
et alB. We found that there was a large disagreement between the two. This comparison is presented in Sec. ||, 



In Sec. II] , we proffer a phenomenological scheme for this underestimation of the band gap and suggest a possible 



reconciliation between the theory and experiment. Conclusions are presented in Sec. LV 



II. COMPARISON : THEORY VS. EXPERIMENT 

In a Letter Ogut et al.a have employed a carefully argued ab-initio methodology to obtain the size-dependent 
optical gap of Si QDs with sizes ranging from 1 to 3 nm.j^These authors jChoose to compare their calculations with 
the experimental observations of Furukawa and MiyasatoQ. A little latcm Buuren et al. reported state of the art 
measurements of the band edges of Si quantum dots (QDs) with average diameters ranging from 1 to 5 nm. By adding 
the measured conduction band (CB) and valence band (VB) shifts to the band gap of bulk Si they obtained the band 
gap of the Si QDs. Recent PL and extended X-ray absorptiori-|data for oxygen terminated silicon nanocrystallites 
of size less than 4 nm, were found to match these observations. These band gaps are smaller than most reported 
theoretical calculations. This factr-was noted by Buuren et al. who choose to compare their observations with older 
calculations by Wang and Zungerll3. We also note that several first principles/LDA calculations on nanocrystalline 
forms of Si have been reported in the pastO. It appears that theoretical calculations over-estimate the gap. 

A comparison of the "state-of-the-art" theoretical calculationsO with the "state-of-the-art" experimental 
observationsQ reveals that the two are in substantial disagreement with each other. This comparison is presented 
in Fig. nl In fact the Ogut et al. calculation does worse than Wang and Zunger (not shown in Fig. m and other 



semi-empirical theoretical calculations. Perhaps the only commonality between the theory and experiment depicted 
in the figure maybe stated in terms of a bland quantum confinement (QC) dictum: "the band gap increases as the 
size decreases" . 

The disagreement between theory and experiment is enhanced at smaller dot sizes. Were we to make a constant 
upward shift of 0.72 eV to the data of Buuren et al. such that they match with the theoretical calculation at 3.5 
nm then the increasing divergence at smaller sizes is clearly manifested. This is depicted in the inset of Fig. |l|. On 
the other hand, we may translate the experimental data horizontally by 1.75 nm and force an agreement with the 
(extrapolated) theoretical calculations. This would imply that the QD sizes have been seriously underestimated by 
Buuren et al.. While the latter do not discount the possibility of some underestimation, a 1.75 nm error is unlikely. 

Ab-initio calculations are computationally demanding at large QD sizes. As we have pointed out the disagreement 
between theory and experiment is pronounced at smaller dot sizes. It should be possible to re-examine or repeat the 
ab-initio calculation in this computationally feasible intermediate regime. 

III. PHENOMENOLOGY 

The observation of visible photoluminescftage-iEL) in a variety of semiconductor nanocrystallites has fueled a large 
body of research work in the past decada3'ElLni3li3. PL has acquired the role of central characterizing tool in this 
field. The photoluminescence spectra from such systems are broad, and often asymmetric about the peak energy. The 
growth of the QDs is a stochastic process. In an earlier work we have argujsd that one needs to consider the distribution 
of crystallite sizes to compare with the experimental spectral lineshapdj'u. A Gaussian size distribution was used, ia 
those works. Fpjt-seHiiconductor nanostructures, the log-normal size distribution has considerable experimentaBil3'E£l 
and theoreticalHila'til support. Specifically, Yorikawa and Muramatsui^ presented an explanation of PL spectra based 
on the log-normal distribution of porous silicon QD size. The experimental work reported by Buuren et al. also shows 
a log-normal distribution. 

We consider the log-normal size distribution, P{d), of diameter d centered around a mean diameter do, 
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dm = exp(^ - a'^) (2) 

do = expfY+^j (3) 

where dm is the dot size for which the maxima occurs in the log- normal distribution and {/i, a} are some characteristic 
constants. 

The number of electrons in a dot of diameter d participating in a PL process is proportional to d^. Thus, for an 
ensemble of QDs, the probability distribution of electrons participating in the PL process is 



Pe{d) = -j= — bd^exp 
\'2TTad 
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where & is a suitable normalization constant. 

In general, the optical band gap is attributed to the energy upshift of the electron and yield 

E = Eo. + §^ (5) 

where E is the enhanced gap, Eoo is the bulk silicon gap (1.17 eV), C is an appropriately dimensioned constant and 
72 is the exponentHu with n e [1,2]. 

Hence the energy upshift, A_E, due to confinement in QD is 

AE^^ (6) 

^Eo = §, (7) 

where AEq is a mean upshift, which is related to the mean diameter do of QD. 



Now, convoluting the upshift (eqn. (^) with the log-normal size distribution (eqn. (||)) 
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The above integral can be easily solved using the property of Dirac delta function, which yields 



d(d) (8) 
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For the PL peak, we equate the derivative of the eqn. (^) to zero. The shift in the PL peak position, A_Ep, is given by 

A£;p = Cexp [-{(3 + n)o-2+^}?i] (10) 



^ / , ^ n[2n+5]/3 



We employ eqn. (ni|) in eqn. (S) to get 



E^Eoo + AEp (12) 
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n[2n+5]/3 



Eq. (13) represents the final form for estimating a more realistic size dependence of the optical gap. The exponenl 
n =1.22 and C = 3.9 in appropriate units, for the results of Ogut et al. shown in fig. (|l|). Kanemitsu et a/.Ej 
obtained the size distribution of oxidized Si nanocrystallites produced by laser breakdown of silane gas. The optimum 
log-normal fit on these experimental data yield dm/ do G [0.7,0.95]. We have used three values 0.7, 0.8 and 0.9 of 
the ratio dm/do in eqn. m3t) and depicted them alongside the data from Ogut et al. and Buuren et al. in fig. (0). 
We see that our calculation with the value 0.7 brings Ogut et aZ.'s result in complete agreement with Buuren et aZ.'s 
result. Further our phenomenological model produces more pronounced downshift at smaller sizes. This is once again 
in agreement with the experimental trend. The asymmetric distribution of the dot size with stretched tailing towards 
larger QDs is responsible for the lower dm/ do ratio, though the skewness may not be as pronounced as dm/da = 0.7. 
Thus, for more asymmetric log-normal distribution one obtains more prominent downshift at smaller sizes. 

IV. CONCLUSION 

We find that a recent "state-of-the-art" theoretical, calculationD for optical gap is in large disagreement with a 
closely following "state-of-the-art" experimental workn. The comparison is even worse for smaller crystallite sizes. 
The experimental work reported a log-normal size distribution for nanocrystallite samples, whereas the theoretical 
work was performed on a single crystallite. This theoretical calculation showed a very good agreement with a decade 
old experimental work, whereas other semi-empirical works were shown to be in poor agreement. To the contrary, we 
find that the semi-empirical works were in better agreement with Buuren et aL's experimental work. 

Ab-initio calculations are computationally demanding at large QD sizes. As we have pointed out the disagreement 
between theory and experiment is pronounced at smaller dot sizes. It should be possible to re-examine or repeat the 
ab-initio calculation in this computationally feasible intermediate regime. 

We have used the methodology discussed in earlier worksaiS with a log-normal distribution. We have shown that 
Ogut et a/.'s results can be in better agreement with experimental works if size distribution is explicitly taken into 
account. In our formulation the downshift of optical gap is larger for smaller crystallite sizes, making this exercise 
even more relevant to the present case. 

We caution however that the disagreement between the theory and the experiment may not be solely related to 
a distribution of crystallite sizes. Other factors could be crucial. These are : (i) partial passivation, (ii) passivation 
with species other than hydrogen, (iii) surface reconstruction, (iv) fiattening of the dot, and-|(v) size underestimation 
by the microscopic techniques. Some of these factors have been mentioned by Buuren et alB. 

There has been a growing realization that the exiatejace. o£,aHpensemble of QDs of varying sizes must be taken into 



account in order to explain experimental observationa3"aE3"IUO'll3. Tentative attempts have been made to understand 
the radiative rate, decay constant, absorption data, etc. based on these considerations. Perhaps any theoretical 
calculation on a single nanocrystallite should be supplemented with effective size averaging before a comparison with 
experiment is made. 
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FIGURE CAPTION 

FIG. 1. The figure compares the "state-of-the-art" experimental band gap obtained by Buuren et al. [Ref. 8] with the 
"state-of-the-art" theoretical one obtained by Ogut et al. [Ref. 6]. The two are in considerable disagreement. In the inset we 
have up-shifted the data from Ref. 8 by 0.72 eV to ensure an agreement with Ref. 6 at 3.5 nm size. The discrepancy between 
the theory and the experiment appears pronounced at smaller sizes. The dashed line is our fit to the calculations of Ref. 6 
{Eg{d) = 1.1 -I- c/d}'"^"^). The '-I-' symbols are data arrived at by using a relationship between CB and VB edges by Buuren et 
al. [see Ref. 8]. 

FIG. 2. The figure depicts the theoretical data of Ref. 6 and the experimental data of Ref. 8. It also shows the band gap of 
crystallites if a log-normal size distribution is incorporated in Ogut et al.'s calculation. Of the four lines depicted in the figure 
the lowest one (dashed line) and the subsequent two (dotted and dashed lines) corresponds to log-normal distributions with : 
dm/do = 0.7, 0.8 and 0.9 respectively. The topmost dotted line is the fit to the data of Ref 6 as mentioned in Fig 1 and the 
text. We notice that the downshift of the band gap is larger for smaller crystallite sizes. We also notice that dm/do = 0.7 
shows excellent agreement with the experimental results of Buuren et al. [Ref. 8] . 
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